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compute  the  kinematic  viscosity  of  water  for  temperatures 
in  the  32-100  degree  Fahrenheit  range  or  0-38  degree  Celsius  range. 

References:  (l)  HDC  Chart  001-1,  Rev  8-6o.  ) 

_ _ —  (2)  Boly,  Ray  E.  and  Tuve,  George  L. ,  ed. ,  CRC  Handbyok  of  Tables 

L,  for  Applied  Engineering  Science,  1970,  page  67^^ 

I  !  !  (3)  Abromowitz,  M.  and  Stegun,  I.  A.,  ed. ,  Handb^0K  of  Mathema¬ 

tical  Functions,  National  Bureau  of  Standa;?<ls  Applied 
Mathematics  Series  55  >  Nov  1970,  page  87'9- 


tmoo  program  is  part  of  the  CORPS  computer  systeiS^ CORPS  is  an 

acronym  for  Conversatioiially  Oriented  Real-Time  Program-Generating  System.  The 
program  is  written  in  time-share  G635  Series,  FORTRAN  IV.  It  consists  of  a 
main  program  and  a  subroutine.  The  main  program  performs  all  I/O  requirements 
with  the  subroutine  performing  all  required  calculations.  All  stops  would  be 
unprogrammed,  resulting  from  improper  input  data.  _ 
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The  program  was  developed  and  is  operational  on  G635,  WES,  Vicksburg,  MS.  It 
is  also  operational  on  HIS  66/80,  Macon,  GA,  and  Boeing  CDC,  Seattle,  WA. 


The  input  variables  will  be  defined  at  execute  time  and  are  the  temperature 
type  (F-Fahrenheit ,  C-Celsius),  number  of  temperatures,  and  the  temperatures. 
The  output  will  consist  of  the  temperatures  as  Fahrenheit  and  Celsius  and 
the  corresponding  kinematic  viscosities  (ft^/sec). 


r.  AOOmOMAL  RIMANKt 

Complete  documentation  of  this  program  is  available  from  the  Engineer 
Computer  Programs  Library,  Technical  Information  Center,  WES. 
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B.  PROGRAM  SPECIFICATIONS; 

Language :  ANSI  FORTRAN  (FORTRAN  IV) 

Solution  Requirements :  The  run  command 

RUN  WESLIB/CORPS/HOOll ,  R 

Method  of  Analysis :  Known  kinematic  viscosities  at  four-degree 
Celsius  intervals  are  stored  in  a  data  statement  and  then  the 
Lagrange  Four  Point  Interpolation  Formiila  is  used  to  determine 
the  viscosities  at  the  input  temperatxires . 

Core  Executive  Requirement:  G635,  10  K  words 

Restrictions :  Temperatures  must  be  in  the  32-100  degree  Fahrenheit 
or  0-38  degree  Celsius  range. 

General  Equations: 

Conversion  of  Celsius  to  Fahrenheit  or  Fahrenheit  to  Celsius 

P  .  |£  .  32  c  . 

where  F  is  degree  Fahrenheit  and  C  is  degree  Celsius 

2 

Kinematic  Viscosity  p  (ft  /sec) 

u  =  _  P(P  -  1)(P  -  2)  ^  (p^.  -..11  f 

^  6-120 

P(P  +  1)(P  -  g)  f  +  P(P.^  1)  r 

2  ^1  6  ^2 

where  f.  ,  for  i  =  -1,  0,  1  and  2  are  known  viscosities  with  f 
being  the  known  viscosity  at  the  nearest  temperature  £  the  given 
temperature  and  p  is  the  ratio  of  the  difference  between  the  given 
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temperature  and  the  temperatiire  at  to  the  difference  between 

the  temperatures  at  f^  and  f^  . 

Range  of  Quantities;  For  temperatures  in  the  32-lC0°F  or  0-38°C 

range;  kinematic  viscosities  will  be  between  1.922  x  10  ^  and 
P 

7.35  X  10“°  ft  /sec. 

Accuracy:  The  kinematic  viscosity  is  computed  to  four  significant 
figures  as  interpolated  from  data  of  the  same  precision  at 
intervals  obtained  from  reference  (2). 
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REF :  ER  1110-1-10-ENGINEERING  AND  DESIGN  -  Engineering  Computer  Prograjn 
Library  Standards  and  Documentation,  Appendix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1 .  PROGRAM  NUMBER:  T22-F3-R0-00K 

2.  TITLE :  HOOll  -  Kinematic  Vicosity  of  Water  for  Temperatures  in  the 

32-100  Degree  Fahrenheit  Range  or  0-38  Degree  Celsius  Range. 

3.  REVISION  LOG:  N/A 

h.  PURPOSE  OF  PROGRAM:  To  compute  the  kinematic  viscosity  of  water  in 

the  32-100°F  or  0-32°C  range. 

References : 

a.  HDC  Chart  001-1,  Rev  8-60. 

b.  Boly,  Ray  E.  and  Tuve,  George  L. ,  ed. ,  CRC  Handbook  of  Tables 
for  Applied  Engineering  Science,  1970,  page  6T. 

c.  Abromowity,  M.  and  Stegun,  I.  A.,  ed. ,  Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards  Applied  Mathematics 
Series  55 >  Nov  1970,  page  879. 

5.  STEP  SOLUTION: 

a.  The  required  inputs,  temperature  type  ( F- Fahrenheit ,  C-Celsius), 
number  of  temperatures,  and  the  temperatures  are  entered  and 
subroutine  HOOll  is  called. 

b.  Temperature  conversions  are  performed. 

c.  The  viscosity  f  at  the  nearest  temperature  ^  the  given  tem¬ 

perature  is  determined  from  the  data  statement  of  known  viscos¬ 
ities.  The  viscosities  ^  i  ^2  then  determined. 

d.  The  ratio  p  ,  the  difference  between  the  given  temperature  and 
the  temperatiire  at  f^  to  the  difference  of  the  temperatures  at 

f,  and  f  is  calculated. 

1  o 

e.  The  kinematic  viscosity  is  determined  using  the  Lagrange  Four 
Point  Interpolation  Formula. 
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f.  The  output,  the  temperatures  as  Fahrenheit  and  Celsius  and  the 
corresponding  kinematic  viscosity,  is  returned  to  the  main 
program  for  appropriate  I/O  operations 

6.  ACCURACY :  The  kinematic  viscosity  is  computed  to  four  significant 

figures  as  interpolated  from  data  of  the  same  precision  at  it°C  intervals 

obtained  from  reference  (b). 


7.  REMARKS: 


a.  Temperatures  must  be  in  the  32-100°F  or  0-38°C  range 

b.  Kinematic  viscosities  will  be  between  1.922  x  10  ^  and 
7.35  X  10“°  ft2/sec. 
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PART  II:  COMPUTER  FUNCTIONAL  DESCRIPTION 


REVISION  LOG:  N/A 
FUNCTIONAL  FLOW  CHART: 
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3.  EQUIPMENT  AND  OPERATING  SYSTEM:  The  program  was  developed  on  a 
G635  time-share  system  in  which  input/output  equipment  consisted  of  a 
Model  33  remote  teletype. 

INPUT  REQUIREMENTS :  All  input  is  entered  via  the  user's  time-share 
terminal  device  in  free  field  format.  The  cues  and  reads  for  input  are 
performed  in  the  main  program.  Since  all  calculations  are  done  in  sub¬ 
routine  HOOll,  all  required  I/O  arguments  for  the  subroutine  are  passed 
via  the  CALL  statement.  The  calling  sequence  for  the  subroutine  is: 

CALL  HOOll ( arg^ » • • . , arg^ , $N ) 

where: 

arg^  -  given  temperature,  degrees  (either  F  or  C) 
arg^  -  type  of  temperature,  F  or  C 
arg^  -  temperature,  °C 

2 

argj^  -  kinematic  viscosity,  ft  /sec 

$N  -  N  is  a  statement  number  in  the  main  program  for  the  non¬ 
standard  return  due  to  temperature  <  0°C  or  >  38°C. 

Arguments  1,  3,  and  1*  are  floating  point,  argument  2  is  character,  and 
N  is  integer.  Arguments  1  and  2  are  the  inputs.  If  argument  1  is 
entered  as  ^C,  it  will  be  returned  as  °F  with  argument  3  as  the  corre¬ 
sponding  °C.  Argument  U  is  the  output  value  of  the  kinematic  viscosity. 
N  is  a  STATEMENT  NUMBER  in  the  CALL  of  the  main  program  which  directs 
return  from  the  subroutine  to  that  STATEMENT  NUMBER  in  the  main  program. 
5.  SECONDARY  STORAGE  INPUT  FORMAT:  None 


U 
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6.  INPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  input 
variables  in  program  HOOll. 

N  -  number  of  temperatures 

temp  -  temperature,  degrees  (either  F  or  C)  (arg^) 

IDEG  -  character,  type  of  temperature,  F  for  °F,  (arg  ) 

C  for  °C 

7.  OUTPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the 
output  variables  in  program  HOOll. 

TEP  -  temperature,  °C  (arg^) 

2 

VISCOS  -  kinematic  viscosity,  ft  /sec  (argj^) 

If  the  input  TEMP  is  entered  as  °C,  it  will  be  returned  from  subroutine 
HOOll  as  °F  with  the  output  TEP  as  the  corresponding  °C. 

8.  PROGRAM  ERROR  MESSAGE:  The  following  is  an  example  of  the  error 
message  printed  from  subroutine  HOOll  if  TEMP  <  0°C  or  >  38°C. 

TEMP  OF  120.50  DEG  NOT  IN  GIVEN  RANGE 

Return  is  to  the  statement  nioraber  N  in  the  main  program  as  described  in 
(!+)  of  this  section. 

9.  VARIABLE  DEFINITIONS: 

HFILE  -  5  character  alphanumeric  program  name;  passed  to  WESLIB 
count  routine  HACCT 

LQZ  -  equal  1,  all  input  cues  and  reads  are  performed;  equal  2, 

call  WESLIB  routine  RERUN  for  entering  only  desired 
input 

LCiX  -  equal  1,  print  instructions  from  RERUN;  equal  3,  no  print 

JKL  -  direct  return  from  RERUN  to  desired  input  read 
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K^.K  -  number  of  input  variables 

IDEG  -  character,  type  of  temperature,  F  for  °F,  C  for  °C 
N  -  number  of  temperatures 

NPOS  -  position  in  DATA  statement  of  the  viscosity  at  the 
nearest  temperature  <  the  input  temperature 

TO  -  temperature  at  NPOS,  °C 

P  -  ratio  of  the  difference  between  the  input  temperature 

and  To  to  the  difference  of  the  temperature  between 
NPOS  +  1  and  TO 


TEMP  -  given  temperature,  degrees  (either  F  or  C) 
TEP  -  temperature,  °C 

2 

V  -  stored  known  kinematic  viscosity,  ft  /sec 

2 

VISCOS  -  kinematic  viscosity,  ft  /sec 


10.  EXAMPLE  CASE:  Determine  the  kinematic  viscosity  of  water  for  6 


given  temperatures. 

a.  Input :  Type  of  temperature  (IDEG)  =  C 
Number  of  temperatures  (N)  =  6 

Temperatures  (TEMP)  =  5»  10,  15 »  20,  25,  30 
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b.  Output : 


INPUT  HOOll  -  KINEMATIC  VISCOSITY  OF  WATER  IN  THE  32-100  DE6-F 
OR  0-38  DEG-C  RANGE 

AA-ENTER  THE  HUMBER  OF  TEMPERATURES  FOR  WHICH  THE  KINEMATIC 
VISCOSITY  IS  DESIRED;HUMBER  MUST  BE  <  OR  =  25 
=  6 

A3-ENTER  TEMPERATURE  TYPE; F=FAHRENHEIT ; C=CELSIUS 
=  C 

AC-ENTER  THE  6  TEMPERATURES  SEPARATED  BY  COMMAS 
=5,10.15,20.5,25.30 


OUTPUT  HOOll  -  KINEMATIC  VISCOSITY  OF  WATER  IN  THE  32-100  DEG-F 
OR  0-38  DEG-C  RANGE 


TEMP 

KIN  VIS 

DEG-F 

DEG-C 

FT>«X2/SEC 

41 .00 

5.00 

0.1634E-04 

50.00 

10.00 

0 . 1407E-04 

59.00 

15.00 

0.1226E-04 

63 .90 

20.50 

0.1067E-04 

77.00 

25.00 

Q.9508E-05 

85.00 

30.00 

0.8518E-05 

RERUN  OR  STOP 
=  STOP 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  III:  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE:  HOOll  -  Kinematic  Viscosity  of  Water  for  Temperatures  in  the 
32-100  Degree  Fahrenheit  Range  or  0-38  Degree  Celsius  Range. 

3.  SOURCE  LISTINGS:  See  pages  9-11. 

U.  NUMERICAL  AND  LOGICAL  ANALYSIS:  Known  kinematic  viscosities  at  four 
degree  Celsius  intervals  are  stored  in  a  data  statement.  The  Lagrange 
Four  Point  Interpolation  Formula  is  used  to  determine  the  viscosities  at 
the  input  temperatures. 

5 .  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  None 

6.  MISCELLANEOUS :  The  program  is  part  of  a  Conversationally  Oriented 
Real-Time  Program-Generating  System  (CORPS).  The  program  is  now  opera¬ 
tional  on  the  WES  0635?  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing 
CDC,  Seattle,  WA.  The  source  listing  on  page  9  contains  the  first  line 
run  command  for  HOOll  and  its  brief.  The  first  line  run  command  runs 
the  binary  HOOllB  of  the  source  listing  on  pages  10  and  11  (the  Fortran 
listing  of  program  HOOll)  and  attaches  the  support  files  RERUN  and 
HACCT. 
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OOOIXSRUN  WE5L IB/CORPS/HO 01  IB, R;WESLIB/RERUN,R:WESLIB/HACCT,R 
0800  5STH1S  PROGRAM  COMPUTES  THE  KINEMATIC  VISCOSITY  OF  WATER  FOR 
0510  59TEMPERATURES  IN  THE  32-100  DEGREE  FAHRENHEIT  OR  0-38  DEGREE 
0820  65CELSIUS  RANGE.  THE  REQUIRED  INPUTS  ARE  THE  TYPE  OF  TEMPERATURE, 
0830  56EITHER  FAHRENHEIT  OR  CELSIUS  AND  THE  TEMPERATURE.  OUTPUT 
OSAO  61C0N5IST  OF  THE  TEMPERATURE  AS  FAHRENHEIT  AND  CELSIUS, PLUS  THE 
0850  3AC0RRESP0NDING  KINEMATIC  VISCOSITY. 

OOSO'-OOPINISH 
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00  0  01)(«RUN  »  =  ;/CORPS/H0  011B(NOGO) 

10000  CHARACTER  HFI L EXS , IDE6»1 
10010  DIMENSION  TP(25) 

10020  HFILE=5HH0011 
10030  LQZ=1;LQX=1 
10040  15000  PRINT  10006 

10050  10006  FORMAT(/''INP:'T  HOOll  -  KINEMATIC  VISCOSITY  OF  WATER  IN  THE  T 
10060*2-100  DEG-F"/16X,''OR  0-38  DEG-C  RANGE"//) 

10080  CALL  HACCT(HFILE) 

10090  GO  T0(15003. 15017), LQZ 
10100  15003  PRINT  15004 

10110  15004  FORMATC'AA-ENTER  THE  NUMBER  OF  TEMPERATURES  FOR  WHICH  THE  KI 
10115*NEMATIC"/"VISCOSITY  IS  DESIRED;NUMBER  MUST  BE  <  OR  =  25") 

10120  15005  READ.N 

10130  GO  10(15007.15017). LQZ 

10140  15007  PRINT  15008 

10150  15008  FORMATC'AB-ENTER  TEMPERATURE  TYPE;F=FAHRENHEIT;C=CELSIUS") 
10160  15009  READ.IDEG 

10170  IF(IDEG.EQ.IHF.OR.IDEG.EQ.IHC)  GO  TO  15012 
10180  PRINT, "ANS  MUST  BE  F  OR  C ; RE-ENTER" ; GO  TO  15009 
10190  15012  GO  T0(15013, 15017), LQZ 
10200  15013  PRINT  15014, N 

10210  15014  FORMATC'AC-ENTER  THE  ",I2,"  TEMPERATURES  SEPARATED  BY  COMMAS 
10215*") 

10220  15015  READ, (TP(I), 1=1, N) 

10230  GO  70(15020,15017), LQZ 

10240  15017  KKK=3 

10250  CALL  RERUN(KKK,LQX, JKL) 

10260  GO  70(15005, 15009, 15015, 15020), JKL 
10270  15020  PRINT  15021 

10280  15021  FORMAT(//"OUTPUT  HOOll  -  KINEMATIC  VISCOSITY  OF  WATER  IN  TH 
10290JE  32-100  DEG-F"/1 7X, "OR  0-38  DEG-C  RANGE"//7X, "TEMP” , 9X, "KIN  VIS" 
103001/"  DEG-F  DEG-C  FTXS2/SEC" ) 

10310  DO  15027  1=1, N 
10320  TEMP=TP(I) 

10330  CALL  H0011(TEMP,IDEG,TEP,VISCOS, $15027) 

10340  PRINT  15028, TEMP, TEP.VISCOS 
10350  15027  CONTINUE 

10360  15028  F0RMAT(F7.2,F9.2,3X,E10.4) 

10370  LQZ=2 

10380  CHARACTER  ZZZZZZ*Z 

10390  16000  PRINT  16002 

10395  16002  FORMAT(/"RERUN  OR  STOP") 

10400  READ  16001,  ZZZZZZ 
10410  16001  FORMAT(A2) 

10420  IF(ZZZZZZ.EQ.2HRE)  GO  TO  15000 
10430  IF(ZZZZZZ.EQ.2HST)  GO  TO  20000 
10440  PRINT, "ERROR  X**  RETYPE" 

10450  GO  TO  16000 
10460  20000  STOP 
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10470  END 

1  1  00  0  SUBROUTINE  HO  0 1 1 ( TEMP , I  DEG, TEP , V I SC05 , *  ) 

11010  CHARACTER  IDEGxi 
11C20  DIMENSION  V(13) 

11030  DATA  V/2.22E-5,1.922E-5.1.686E-5,1.491E-5,1.3JE-5,1.194E-5, 
110404  1 .08E-5,9.829E-6.8.993E-6,8.269E-6,7.636E-6,7.08E-6,6.589E-6/ 
11050  IFdDEG.EQ.  IHC)  TEP  =  TEMP 
11G60  IF(IDEG.EQ.IHC)  TEMP=9 . *TEP/5 . +32 . 

1  1  070  IF(IDEG.EQ.IHF)  T  EP  =  (  T  EMP-32  .  )  5  . /9  . 

IIOCO  IFCTEP.GE.O. .AND.TEP.LE.38.)  GO  TO  10110 
1  1  0  9  0  IF( IDEG. EQ . IHF)  PRINT  1 0 0 90  , TEMP , I  DEG 
11100  IFCIDEG.EQ. IHC)  PRINT  10090 , TEP, IDEG 

1  1  1  10  1  00  90  FORMATC'TEMP  OF  ".F6.2,"  DEG-",A1.’'  NOT  IN  GIVEN  RANGE") 

11120  RETURN  1 

11130  10110  NP0S=TEP/4+2 

11140  TO=NPOSX4 . -8 . 

11150  IF(TO-TEP)  10150,10140,10140 
11160  10140  VISC03=V(NP0S);RETURN 
11170  10150  P:(TEP-T0)/4. 

11120  VI5COS:V(NPOS-l  )  ( P-P*P  ) )♦  ( P-2  .  )/6  .  +  V {  NPOS )  «  (  PxP-1  .  ) 
11190*x(P-2.)/2.  +  V(NP0S  +  1  )x<2.  «P-P!<P)»«(P+1  .  )/2  . 

]  120Ci  +  V(NPOS  +  2)»<(P>  P<<P-P)/6  . 

11210  RETURN 
11220  END 


li 


